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Abstract 

We investigate p-wave pion production in nucleon-nucleon collisions up to next-to-next-to-leading 
order in chiral effective field theory. In particular, we show that it is possible to describe simulta- 
neously the p-wave amplitudes in the pn — ► ppir~, pp — > pmr + , pp — > dir + channels by adjusting a 
single low-energy constant accompanying the short-range operator which is available at this order. 
This study provides a non-trivial test of the applicability of chiral effective field theory to reactions 
of the type NN -> iViVvr. 

PACS numbers: 11.30.Rd, 12.39.Fe, 13.60.Le, 21.30.Fe, 25.10,+s 
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I. INTRODUCTION 



With the advent of chiral perturbation theory (ChPT), the low-energy effective field 
theory (EFT) of QCD, high accuracy calculations for hadronic reactions with a controlled 
error estimation have become possible [H, 0] . In that framework, tttt [3| and ttN jij scattering 
observables and nuclear forces [§J are calculated based on a perturbative expansion in q/A x 
with q referring to either a generic momentum of external particles or the pion mass m n , 
and A x ~ 1 GeV being the chiral symmetry breaking scale. An extension of this scheme 
to pion production in nucleon-nucleon (AW) collisions turned out to be considerably more 
difficult. A straightforward application of the power counting proposed by Weinberg 0, 0] 
to the reactions AW — > NNtt @, Q failed badly (see also Ref. [lOj where it was pointed 
out that the naive power counting using the heavy baryon formalism is not applicable above 
the pion production threshold). Indeed, for neutral pion production in pp collisions, the 
corrections due to the next-to-leading order (NLO) increased the discrepancy with the data 
and, moreover, the next-to-next-to-leading order (NNLO) contributions turned out to be 
even larger than the NLO terms [11]. The origin of these difficulties was identified quite 
early by Cohen et al. see also jl3j]. who stressed that the additional new scale, inherent 
in reactions of the type AW — > NNtt, needs to be accounted for in the power counting. 
Since the two nucleons in the initial state need to have sufficiently high kinetic energy to 
produce the onshell pion in the final state, the initial center-of-mass momentum needs to be 
larger than 

"' (1) 



Pthi = \jM N mi< 



with 



Pthr 

A. 



0.4, 



where m n and Mn refer to the the pion and nucleon mass, respectively. The proper way to 
include this scale was presented in Ref. Q and implemented in Ref. [l5j], see Ref. [l6[ for 
a review article. As a result, pion p-wave production is governed by the tree- level diagrams 
up to NNLO in the modified power counting scheme of Ref. [13]. On the other hand, for 
pion s-wave production, pion loops start to contribute already at NLO. It was demonstrated 
in Ref. TtJ that all irreducible loop contributions at NLO cancel altogether, and the net 
effect of going to NLO was shown to increase the most important operator for charged pion 
production, first investigated in Ref. fl~8l |. by a factor of 4/3. This was sufficient to overcome 
the apparent discrepancy with the data in that channel. But the neutral pion channel is 
more challenging — it still calls for a calculation of subleading loop contributions. First 
steps in this direction were taken in Refs. We further emphasize that the A(1232) 

isobar should be taken into account explicitly as a dynamical degree of freedom 12| because 
the Delta-nucleon mass difference, AM, is also of the order of pth. r - This general argument 



21, 22 



was confirmed numerically in phenomenological calculations [2( 

Pion p-wave production in AW collisions receives an important contribution from the 
leading (NN) 2 tt contact term in the effective Lagrangian, which also figures importantly 
in the three-nucleon force 
processes 'yd — *> ttNN [24j 
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231 ] . In addition, the same operator also contributes to the 

as well as to weak reactions such 
21, 
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as visualized in Fig. [TJ 



and ird -> 7 AW j; 
as, e.g., tritium beta decay and proton-proton (pp) fusion [2i 
Notice that this operator appears in the above reactions in very different kinematics, ranging 
from very low energies for both incoming and outgoing AW pairs in pd scattering and the 
weak reactions up to relatively high initial energies for the AW induced pion production. 



In Ref. 30j it was shown that both the 3 H and 3 He binding energies and the triton /3-decay 



can be described with the same contact term. However, an apparent discrepancy between 
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FIG. 1: Illustration of the various reactions, where the leading (iViV) 2 7r-contact term contributes. 

the strength of the contact term needed in pp — > pmr + and in pp — > de + v e was reported in 
Ref. [31]. If the latter observation were true, it would certainly question the applicability of 
chiral EFT to the reactions NN -> iViVTr. 

To better understand the discrepancy reported in Ref. |3lj, in this paper we simulta- 
neously analyze different pion production channels. In particular, we calculate the p-wave 
amplitudes for the reactions pn — > pprc~ , pp — > pmr + , and pp — > d-K + . Note that even in 
these channels the contact term occurs in entirely different dynamical regimes. For the first 
channel p-wave pion production goes along with the slowly moving protons in the x Sq final 
state whereas for the other two channels the pp state is to be evaluated at the relatively 
large initial momentum. Notwithstanding, all three channels of the reaction NN —>■ NNn 
seem to give consistent results for the low-energy constant (LEC) d that represents the 
strength of the contact term, as we will show in the present paper. We discuss which ad- 
ditional data are needed to further support this conclusion. We argue that the origin of 
the discrepancy reported in Ref. 131] is not due to the different kinematics of NN — > NNir 



and pp fusion, but rather in the inconsistency in the partial wave amplitudes used in the 



analysis. In addition, we also comment on technical issues related to the work of Ref. [31 

Our manuscript is organized as follows: In Sec. [Ill we discuss the general features and the 
relevant observables for p-wave pion production. In Sec. [Till the power counting is outlined 
with special emphasis on the p-wave amplitudes. Our results for the various pion production 
channels are presented in Sec. IIV1 Here, we also discuss the role of the leading ttN scattering 
parameters, C3 and C4, for the p-wave pion production amplitudes. We close with a short 
Summary. 



II. GENERAL REMARKS 

It is not obvious, a priori, that with just a single contact term, which contributes to 
the various reactions shown in Fig. [U a consistent description of all these channels can 
be achieved. The purpose of the contact term is twofold: it should, on the one hand, 
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absorb any sensitivities to the employed NN wave functions and in this way remove the 
model dependence in the evaluation of the observables. On the other hand, it provides a 
parameterization of the short-range physics that contributes to the process being considered. 
Thus, the strength of the contact term is necessarily dependent on the method applied to 
regularize the integrals (typically a cut-off) and also on the NN interaction that is used for 
generating the wave functions. 

For the case at hand the contact term connects NN S'-waves in the initial state with 
NN S-waves in the final state. Since the contact term is a local four-nucleon operator, after 
including the NN distortions its contribution scales as the product of the initial and final 
NN wave functions at the origin. Each of these wave functions, in turn, may be represented 



by the inverse of the corresponding Jost function [32|. The reason why it is expected to 
be possible that the same contact term can be used in all reactions listed above is that the 
energy dependence of the Jost function is fixed by the onshell NN data and is therefore 
independent of the unknown short range physics. Specifically, the NN distortions can be 
represented as an integral over the relevant phase shifts by means of the so-called Omnes 
function [33| — see also the discussion in Ref . [13] • This is correct up to contributions from 
the left-hand cuts and the high energy behaviour of the NN interaction, both are expected 
to be of higher order in the expansion. As opposed to the energy dependence, the overall 
scale of the distortions can be shown to be sensitive to things like the NN interaction and 
the cut-off employed j35j. Clearly, what needs to be assumed in the argument given is that 
there is a proper separation of scales in the problem. Note that the expansion parameter 
p t hr/A x ~ 0.4 is quite large. In this sense a consistent description of all mentioned reactions 
with the same contact term provides a non-trivial test of the applicability of the chiral 
expansion to pion production in NN collisions. 

One might ask why we take the effort of this study, since in Ref. [3l[ it was already shown 
that a consistent description is not possible. The answer is twofold: first of all, we found 



that the partial wave decomposition of Ref. [36| , the result of which was used in Ref. (31 
is not correct (see discussion in Sec II V C|) — this is why we decided to directly compare to 
the data in the present work. Secondly, there is also a conceptual problem in the work of 



Ref. [31|: as was outlined above, as long as different phase-equivalent NN interactions are 
used, it should be possible to absorb the model dependence of the calculation in a single 
counter term up to higher-order corrections. However, in Ref. [31] pion production from 
initial NN and NA states is not treated on equal footing. Rather the contribution from the 
A isobar excitation is added on top of and independently of the employed NN interactions. 
Thus, it is quite possible that the utilized NN — > NA transition potential is too strong. 
Specifically, it is not constrained by the empirical NN phase shifts as it is the case when 
considering the NN and NA amplitudes consistently within a coupled-channel (Lippmann- 
Schwinger-like) scattering equation [2(j. In this sense, it should not come as a surprise that 
it was not possible to absorb the model dependencies in a single counter term within the 
scheme used in Ref. |3llj. To avoid this problem, in this work we employ the coupled-channel 
NN model of Ref. |37| which involves the NN — > iVA transition potential. 

Eventually all reactions shown in Fig. [1] should be analysed consistently. This would, 
however, require a calculation to third order in the chiral expansion of the process 7<i — > 
ttNN and rrd — > '-/NN or a rather involved three-nucleon calculation for the tritium beta 
decay which goes beyond the scope of this work. Instead, as a next step in this ambitious 
program, we analyse here in detail various pion production channels. Notice that although 
these reactions appear to have similar kinematics, the relevant transition for the reaction 



4 



pn — > ppir~~ involves very low momenta in the NN ^So state and considerably higher momenta 
(~ Pthr) i n the 3 Si channel while the situation is just opposite for the reactions pp —* 
(d/pn)iT + . Thus, a simultaneous description of these reaction with a single short-range 
operator indeed provides a highly nontrivial consistency test of our approach. Notice further 
that the (NN) 2 tt short-range operator we are interested in here does not contribute to the 
pp —>■ pp7t° transition which is, therefore, not considered in the present work. Thus, the only 
reactions of interest for this study are pp — > (pn/d)n + and pn — > ppn~ . Here the relevant 
transitions are 1 Sq — ► 3 S±p for the former and ( 3 S\ — 3 Di) — > ^SqP for the latter, where the 
small letter labels the pion angular momentum. Since the main focus of this work is on the 
role of the contact term, we will concentrate on observables where the final NN system is in 
an S-wave — which largely simplifies the numerical work. However, as outlined below, the 
contribution of NN P- waves to observables might be relevant for the reaction pp — > pmi + . 
This potential problem renders this channel not very convenient for the extraction of the 
counter term, as will be discussed in Sec. HVl 

To be specific, we calculate in this work the differential cross sections and analyzing 
powers for the reactions pp — > d7c + , pp — > pnn + , and pn — > (pp)is () 7i~ ■ Here the symbol 
(pp)is indicates that in the corresponding measurement the final pp relative momentum 
was restricted kinematically to be less than 38 MeV/c (M pp — 2M^ < 1.5 MeV) which leads 
to a projection on the x So pp final state. For all the mentioned observables experimental 
data are available or will be available soon in the energy range of relevance here. Besides the 
anisotropy of the pion angular distributions, all observables are sensitive to both s- and p- 
wave pion production. Although there exists an NLO calculation for s-wave pion production 



in pp — > dn + using ChPT, its theoretical uncertainty is still sizable [12J. For s-wave pion 
production accompanied by a transition of an isospin-one NN pair to an isospin-one NN 
pair (e.g. in pp — > ppir ), no sufficiently accurate ChPT calculation is available at present. 
Since we focus here on the p-wave amplitudes, we extract the s-wave amplitudes directly 
from the data in order to minimize the uncertainties of our calculation. The phase of these 



amplitudes is then imposed using the Watson theorem [32|, see the discussion in Sec. HVl 



It is well known that p-wave pion production in pp — > dn + and in pp — > pmr + is strongly 
dominated by the transition 1 D 2 — ► 3 Sip due to a strong coupling of the initial iViV state 



to the 5 S2 NA state 38j. Therefore, the amplitude we are interested in has only a minor 
impact on the observables. In other words, the uncertainty for the extraction of the counter 
term from these reactions will be significant. The situation is much more promising for the 
reaction pn ppir~: here the amplitude of interest is the leading p-wave. In addition, the 
strength of the s-wave amplitude can be taken from the reaction pp ppir using isospin 
symmetry and correcting for the final state interaction (FSI) as discussed in Sec. IIVDI 
Unfortunately, no data are presently available for pn — > ppir~ at sufficiently low energies. 
Nevertheless, as will be shown below, already the higher-energy data provide some insights. 
In addition, data at lower excess energies will be available soon [39]. 

The goal of the present investigation is to explore whether it is possible to obtain a 
simultaneous description of all iViV — > NNtc channels. A more quantitative study including 
a statistical analysis of the data and an estimation of the theoretical uncertainty is postponed 
until accurate experimental data will become available for pn —>■ ppir~ at low energies. 
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III. FORMALISM 



Our calculations are based on the effective chiral Lagrangian with explicit A degrees of 



freedom. The leading irN and 7riVA interaction terms read |40l . 141 



£(0) = N i 



1 



9 A ^ 



h A r 
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+■ ■ 



(2) 



while the first corrections have the form 
1 



iN ] r ■ (tt x Vtt) • VA + /i.e.) - —A* 



c 3 (V7r) 2 + - c 4 + 



4M 
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A - t N\t ■ a ■ Vtt)A A f A + 



(3) 



The ellipses stand for further terms which are not relevant for the present study. In the 
equations above f n denotes the pion decay constant in the chiral limit, g& is the axial- vector 
coupling of the nucleon, Iia is the AA-7T coupling, iV and \I/a correspond to the nucleon 
and Delta fields, respectively, and S and T are the transition spin and isospin matrices, 
normalized according to: 



S,S] 



~(25ij - ie ijk T k ) 



(4) 



We also emphasize that the effective Lagrangian of Ref. 40] contains another (NN) 2 tt con- 
tact operator which can be shown to be redundant as a consequence of the Pauli principle 



u, m 
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We are now in the position to discuss the relevant scales and counting rules for p-wave pion 
production. We assign the outgoing two-nucleon relative momentum p' and the outgoing 
pion momentum k n to be of order of and introduce the expansion parameter 
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AM 
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(5) 



where p ~ p^ is the initial two-nucleon relative momentum. The counting rules for the time- 
dependent vertices, such as e.g. the Weinberg-Tomosawa (WT) vertex in are discussed 
in detail in Refs. [13, At leading order one finds that the WT vertex is oc 2u) m with u n 
being the energy of the outgoing (onshell) pion. The diagrams contributing to the production 
operator at LO and at NLO are shown in Fig. [2] whereas the corresponding graphs at NNLO 
are depicted in Fig. [3j At NLO there are only diagrams in which the pion is produced 
through the excitation of the A resonance. The relative suppression of these diagrams as 
compared to the ones involving the nucleon is accounted for by the A propagator which is 
suppressed by 1/p as compared to l/m n in the nucleon case. To see that the diagrams in 
Fig. [3] indeed contribute at NNLO for p-wave pion production consider, as an example, the 
first graph in this figure. Its contribution can be estimated using dimensional analysis as 
follows: 



(6) 
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FIG. 2: " (Color online)" Leading and next-to-leading diagrams for the p-wave amplitudes of NN — > 
NNtt. Single (double) solid lines denote nucleons (Deltas), dashed lines denote pions, green ellipses 
correspond to the NN wave functions in the initial and final states. 




NLO 



FIG. 3: "(Color online)" Diagrams that contribute at NNLO to the p-wave amplitudes of NN 
NNtt. Subleading vertices are marked as 0. 



Here we used that the outgoing pion momentum k w enters the ttNN vertex to allow for the 
p-wave amplitude. To understand the suppression factor this operator should be compared 
with the LO contribution k^ / (f^m^) . Thus, one gets an order x 2 suppression for the first 
diagram of Fig. [31 Similarly, using the ttttNN vertex from £W i n combination with the p/ f n - 
scaling for the ttNN vertex one arrives again at a x 2 suppression for the second diagram in 
Fig. [3J Further details can be found in Appendix [X] which contains explicit expressions for 
the diagrams shown in Figs. [2] and El Once the amplitudes are evaluated they need to be 
convoluted with proper NN wave functions. Ideally, one would use wave functions derived 
from the same formalism, namely ChPT. However, up to now these are only available for 
energies below the pion production threshold We therefore use the so-called hybrid 
approach, first introduced by Weinberg [7J, based on the transition operators derived within 
the effective field theory and convoluted with realistic wave functions 37J]. This procedure 
should also provide reasonable results, however, a reliable uncertainty estimate is possible 
only at the level of the transition operator. 
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IV. RESULTS AND DISCUSSION 



A. Parameters of the calculation 

To the order we are working, the following low-energy constants (LECs) appear in the 
calculation: f w , g A , h A , C3, C4 and d. Only the last LEC cannot be taken from other sources, 
for its value strongly depends on the NN wave functions employed. We adopt the following 
values of the parameters: f n = 92.4 MeV, g A = 1.32, h A ~ 2.1g A = 2.77, c 3 = -0.79 GeV^ 1 
and C4 = 1.33 GeV^ 1 . The values of the LECs C3 and C4 are taken from Ref. j43|. From the 



fit to iiN threshold parameters, two solutions for the q are given in Ref. [43| corresponding 
to the different choices of h A {h A ~ 2.1g A and h A ~ 2.1). The sensitivity of the results 
to the different values of C3 and C4 will be also discussed. As already mentioned in the 
Introduction, the power counting scheme calls for a dynamical treatment of the A isobar 
as a result of the comparable numerical value of the Delta-nucleon mass difference and pthr- 
The implications of integrating out the A degrees of freedom for the processes at hand are 
discussed in Appendix iBl 

The deuteron wave function and the NN scattering amplitudes used in the calculation are 
generated from the CCF NN potential [33] - As described above, we do not calculate the s- 
wave pion amplitudes in this work but rather take both their strength and the phases directly 
from experiment. To be specific, for the reaction pn — > ppn~ we aim at the description of the 



double differential cross sections and the analyzing power measured at TRIUMF [44], |45j and 



PSI 1461. Following the Watson theorem to parameterize the relevant 3 Pq — > ^Sqs amplitude, 
we use the ansatz Ce lSsp o^^ (r = 0), where the inverse Jost function in the ^So partial 



wave, (r = 0), and the initial phase shift 5i Po are calculated from the NN model used, 
and the parameter C is fitted to reproduce the corresponding amplitude extracted from 



'esp 

the TRIUMF data using a partial wave analysis [45|, 47j j. It is interesting to note that the 
3 P -> l S Q s amplitude from the TRIUMF analysis at 1 r/=0.66 (T lab = 353 MeVHs about 25% 
larger than that extracted from the pp — > ppn° measurement at CELSIUS 48J. A similar 



inconsistency is discussed in Ref. |49j] where it is argued that the total cross sections at low 
energies for pp — > ppir recently measured at COSY are about 50% larger than those at 
CELSIUS and IUCF as a result of the missing acceptance at small angles for both CELSIUS 
and IUCF. 

For the reaction pp — > d-K + the s-wave amplitude occurs in the 3 P\ — > 3 Sis partial wave 
and can be related to the total cross section at threshold. The most precise way of getting 
this quantity is to extract it from the width of pionic deuterium atom, measured at PSI 



with high accuracy [50j, ]5JJ. This procedure gives the following value of a, the total cross 
section divided by r\: a = 252^ fib j52|. Thus, we adjust the mag nitude of the 3 Pi -> 3 Sis 
amplitude to be in agreement with this observable. 

As mentioned above the value of d depends on the iViV interaction employed and on the 



method used to regularize the overlap integrals. Indeed, in Refs. [141 . |28| a strong sensitivity 
of the LEC d to the regulator is reported. It therefore does not make much sense to compare 
values for d as found in different calculations. What makes sense, however, is to compare 
results on the level of observables and this is what we will do below. We will adjust the 



1 Traditionally, the energy in the pion production reactions is given in terms of 77, the (maximum) pion 
momentum allowed in units of the pion mass. 
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FIG. 4: "(Color online)" Results for A 2 /Aq (see Eq. [7J (left panel) and the analyzing power at 
90 degrees (right panel) for the reaction pp — ► dir + for different values for the strength of d (in 
units 1/(7^ Mat)). Shown are d = 3 (red solid line), d = (black dashed line), and d = —3 (blue 
dot-dashed line). The data are from Refs. [53j, [54J, |55j, 
s-wave amplitude is fixed from data. 



56, 



571 ]. The strength and phase of the 



value of d in such a way to get the best simultaneous qualitative description of all channels 
of NNn. 



B. Reaction pp — > d-K + 

We begin with a discussion of the results for the reaction pp — > dix + . In Fig. HJ we 
compare our calculation for various values of d with the experimentally available angular 
asymmetry parameter A 2 /A . The coefficients A t are related to the unpolarized differential 
cross section via 

^ = A Q + A 2 P 2 (cos9 7T ) , (7) 

with P 2 (x) being the second Legendre polynomial and 6 n the pion angle in the cm. frame. 
We also show results for the analyzing power at 90 degrees. In both cases the observables 
are plotted as functions of the parameter 77. Here and in what follows, the value of the LEC 
d is always given in units l/(/^Mjy). Notice that at low energies, it is sufficient to just show 
the analyzing power at 90 degrees since its angular dependence is proportional to sin 6^. To 
illustrate this we also present in Fig. [5] the analyzing power as a function of the scattering 
angle for two different energies 77=0.14 and 77=0.21. At rj ~ 0.5 the angular dependence of 
the analyzing power starts to deviate significantly from sin 8 n due to the onset of (i-waves. 
Clearly, at these (and higher) energies we cannot expect our calculation to agree with the 
data anymore. As can be seen from the figures, the data at small rj, especially the analyzing 
power, prefers a positive value for d — our fit resulted in d — 3 for the best value. To 
demonstrate the effect of the LEC d on the observables, in Fig. and in subsequent figures 
we also give the results with d = and with the negative LEC d =-3. 
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180 



FIG. 5: "(Color online)" Re- 
sults for the analyzing power 
at 77=0.14 (upper panel) and 
r/=0.21 (lower panel) as func- 
tions of the angle 9 n for the re- 
action pp — * dir + for different 
values of d. Shown are d = 3 
(red solid line), d = (black 
dashed line), and d = —3 
(blue dot-dashed line). The 



data are from Ref. 56f] . The 
strength and phase of the s- 
wave amplitude is fixed from 
data. 



e, (cleg) 



C. Reaction pn — ► ppit 

We now turn to the reaction pn — > ppir~ . As it was already explained in Sec. [Til f° r 
this reaction channel the relevant pion p-wave occurs in conjunction with the two-nucleon 
pair in the ^So state. It is known experimentally that in the isospin-one channel the final 
P-wave diproton contributions (Pp and Ps) start growing with the energy rather rapidly 
so that already for excess energies around 30 MeV they provide about 50% of the total 
cross section [48|. Therefore, in order to be sensitive to our particular amplitude one needs 
to isolate experimentally the S'-wave diproton state by putting kinematical cuts on the 
two-nucleon relative momentum. This is exactly what was done in the experimental study 
of pn — > ppir~ at TRIUMF 0, H. In particular, they measured the differential cross 



section d a/(dVtdm pp ) and analyzing power A y for T iab = 353 MeV (77=0.66), where the final 



diproton relative momentum p' was restricted to be not larger than 38 MeV/c (M pp — 2Mm ~ 
1.5 MeV). A similar measurement for the analyzing power was also performed at PSI |46j for 
Tiab = 345 MeV and pp invariant masses M pp — 2M N < 6 MeV. It is interesting to note that 
the positions of the peaks in A y seem to be somewhat different in these experiments (se e 
Fig- E]), although the data of Ref. 46[ have much larger uncertainties than those of Ref. 44 



Unfortunately, presently data for pn — > ppn~ are only available at such high energies where 
our corresponding results in the pp — > d-K + channel already start to deviate considerably from 
the experiment. Therefore, for the reaction pn — > ppix~ we expect likewise only a qualitative 
description. Nevertheless, a comparison with the experimental data in this channel is quite 
instructive too and shows also a preference for a positive value of d as visualized in Fig. 
Fortunately, there will soon be a measurement for the same observables at lower energies at 
COSY [3^]. Once this data will be available we should be able to draw more quantitative 
conclusions on the value of the parameter d needed for the reaction pn — ► ppn~ . 



10 




60 120 

0„ (deg) 



60 120 

Q (deg) 



180 



FIG. 6: "(Color online)" Results for d 2 a/d^l n dMp p (left panel) and A y (right panel) for pn — ► 
pp( 1 5o)vr~. Shown are the results for d = 3 (red solid line), d = (black dashed line) and d = —3 



(blue dot-dashed line). The data is from TRIUMF 
circles) . 



44 



451 ] (black squares) and from PSI [46j] (blue 



D. Reaction pp — ► pmr 



The reaction pp — > pmr + is the most difficult and the least convenient one for the ex- 
traction of the contact term. Besides the fact that here, as in pp — ► dir + , pion p-wave 
production is mainly driven by the l D 2 initial state, in addition NN P-waves contribute 
for isospin-one as well as for isospin-zero NN final states. At the energies considered in 
the experimental investigation, rj =0.22, 0.42, and 0.5, the Pp amplitudes may contribute 



significantly |48|, |58|, 159 



They should be particularly important in view of the smallness 
of the amplitude — even small contributions to A 2 , see Eqs. dZj) and can affect the 
partial wave analysis considerably. In the partial wave analysis performed in Ref . 36[ , these 
Pp contributions were not taken into account at all. Also there are contributions to Aq from 



the isospin-one NN final states that potentially increase the uncertainty of the analysis, 
especially in view of the differences in the experimental results in pp — >• ppn as already 
discussed above. These arguments alone cast serious concerns on the partial wave analysis 
performed in Ref. [36]. But there is an even more direct evidence of problems with the 
extraction of the partial wave amplitudes of Ref. 36[ which we now discuss in detail. The 
observables measured for the reaction pp — - > pmc + in Ref. [36] include the coefficients A and 
A 2 in the differential cross section, see Eq. (j7]) and the analyzing power Ay (90°). Neglecting 
the Pp contributions these observables can be expressed in terms of the three partial wave 
amplitudes with the isospin-zero pn-state ao (^o — > 3 Sip) (the single amplitude, where the 
(NN) 2 tt contact term contributes), cq ( 3 Pi — > 3 5*is), a 2 (^D 2 3 Sip) and the contribution 
of the isospin-one channel denoted as A 1 ^ 1 via 



A y (90°)[A 



A 

A 2 
A 2 



Idol 2 + ki| 2 + I CX 2 1 2 



4 

-(v^ImfaiOo] + ImfaiG^]), 



(8) 
(9) 
(10) 
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FIG. 7: "(Color online)" Results for the magnitude of A 2 (left panel) and A y (90°)(A - A 2 /2) 
(right panel) for the reaction pp — > pnir + for different values of the contact term. The notation of 
curves is the same as in Fig. HJ The data are from Ref . 36 ] . 



where (A - A 2 /2) is just da/dtt(90°) from Eq. Using the system of Eqs. fjH])- (flT)]) one 
can determine the amplitudes ao, ai and a 2 provided one knows the isospin-one piece Aq =1 . 
The latter was extracted in Ref. 13a] from the measurement of the total cross section in the 



reaction pp — ► ppn reported in Ref. 60]. However, the FSI in the pp —>■ ppn^ reaction is 



o 



very different to that in the pp —>■ pn7c + channel. To estimate the difference note that in 
the energy region studied, which is less than 20 MeV, the dominant partial wave is the one 
where the final two-nucleon state is in the S'-wave. In this case the correction factor would 
be proportional to the ratio of the inverse Jost functions squared integrated over the phase 
space 

Jd 3 p'k w \F pn (p')\ 2 

where F pn (p') and F^^p') are the inverse Jost functions for the pn and pp 1 So states, respec- 
tively. As discussed above, although the Jost function itself depends on the NN model used, 
its energy dependence does not. We may, therefore, evaluate R using any sensible model 
for the NN interaction. For a separable NN potential there exists an analytic ex pre ssion 



for the Jost function in the pp system in the presence of the Coulomb interaction [6l|, [62 



Using it one finds that the ratio R is about 1.5 for rj = 0.22 and about 1.2 for rj = 0.42 



Similar results are obtained using the CCF NN interaction [37J. Thus, compared to the 



original analysis performed in Ref. [36J], the isospin-one contribution at r\ = 0.22 should be 
enhanced by more than a factor of two if, in addition, one utilizes the new, larger experimen- 
tal data from COSY for the total cross section for pp — > ppw |4{|. This change, of course, 



will significantly affect the results of the partial wave analysis. Given the above difficulties 



with the partial wave analysis of Ref. [36j, we decided to compare our results directly to the 
experimentally measured quantities. Aiming presently at a qualitative description of the 
data, we will not include the Pp-states in this work. 

The results of our calculation for A 2 are shown in Fig. [3 Again, positive values of the 
contact term with d ~ 3 seem to be preferred. We emphasize, however, that these results 
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0.2 0.3 0.4 0.5 



FIG. 8: "(Color online)" 
Comparison of ao as it results 
from our analysis in compar- 
ison to the partial wave am- 
plitude extracted in Ref. [36f] . 
The notation of curves is the 
same as in Fig. |U Since our 
calculations well describe all 
observables of Ref. [30], this 



figure nicely illustrates the 
problem of the partial wave 
decomposition of this refer- 



ence. 



should be treated with care since the calculations at higher energies can be affected by P- 
wave contributions whereas the lowest point is not very sensitive to the value of d due to the 
large experimental uncertainty. We can also check whether our results are consistent with the 
measurement of the analyzing power that is related to our amplitudes via Eq. ({TO]) . To allow 
for this comparison, however, we need to know the pion s-wave amplitude a\. At present, 
this quantity is known theoretically only up-to-and-including terms at NLO. Therefore, to 
minimize the uncertainty of the current study, we extract this amplitude directly from data 
on the total cross section in pp — > pnir + through Eq. (JSJ. We employ the amplitude Aq =1 
consistent with the data at COSY and correct for the FSI factor as described above and 
take the amplitudes ao and a 2 from our NNLO calculation. In the right panel of Fig. [7J 
we compare our results for A y (90°) (A — A 2 /2) with the corresponding data. Since we use 
the experimental total cross sections to extract ai, our results, given by red (d — 3), black 
(d = 0) and blue (d = —3) squares in the right panel of Fig. [JJ can be presented at specific 
energies only. The squares include the experimental uncertainty in the total cross section 
A Q used to extract a\. To guide the eye, we also show the results of interpolations between 
the three energies. It is seen that the magnitude v4 y (90°) (A — A 2 /2) is much less sensitive 
to the value of the LEC d than, e.g., A 2 . Notice further that the experimental points do not 
include a 12% uncertainty due to systematic errors in A and A 2 . 

Since we do not know the contribution of the NN P-waves to the pp — > pmr + observables 
at present, and an improved partial wave analysis would require a careful study of various 
uncertainties, we do not try to extract ao from the data. However, in order to illustrate the 
potential effect of the changes discussed above (up to NN P-waves) on ao, in Fig. [8] we show 



the results of our calculation for ao in comparison to the old extraction of Ref. 36]. Evidently, 
although all data presented in Ref. [36| are in a good agreement with our calculation (as 
demonstrated in Fig. [7]), the partial wave amplitude is not at all described — see solid curve 
for our results with d=3 in Fig. [HI which illustrates clearly that the partial wave solution given 



in Ref. [36] should be abandoned. It is interesting to note that in Ref. [31] it was stressed 
that a positive value for ao is necessary in order to achieve a result for pion production that 
is consistent with the ones for the weak rates. This is in accord with our findings based 
solely on the data for NN — > NNtt. Here we do not aim at a more quantitative comparison 
with Ref. jUJ, because of the technicalities discussed in the beginning of Sec. [UJ 



Finally, we would like to discuss the sensitivity of our results to the parameters q. As 
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shown in Appendix [S] for the 1 S , — > 3 Si or 3 Si — > ^So iViV transitions the parameters q 
occur in the combination Cf 1 = c%/2 + C4 + l/(4Mjy). This combination appears to be 
largely constrained by the ttN data since the different sets of q from the recent analysis [43[ 

3c- 

give basically the same value for C i 1 . In addition, to the order we are working at, this 
combination is fully absorbed in the counter term since the corresponding potential for 
NN — > NNir, see the second diagram in Fig. [31 is just a constant up to higher order terms 

~ ^ V^rT"? = C? '' ( l + 0(x2) ) ' (12) 

Due to a coupled channel effect, the same combination of q also contributes in the 3 Di — > 
3 Si — ► partial wave. The situation is different when D-waves contribute at the level of 
the transition operator. The combinations of q in the 3 Di — > l SoP amplitude for pn — > pp7r~ 
and in the 1 D 2 — > 3 Sip amplitude for pp —>■ (d/pn)n + will influence the observables, for at the 
order we are working at there is no contact term that can absorb the resulting dependence 
on the LECs q. It is worth mentioning that the combinations of q in these partial waves are 
constrained only weakly by nN data. In particular, the combination of q in the 3 D\ partial 
wave, C^ Dl — c.s — c 4 — 1/(4M N ), changes from 2 to 7 depending on which of the sets of c$ 
given in Ref. [43| is used. Thus, we conclude that the reaction AW — > NNn may serve as an 



additional source of information to constrain the q's, complementary to the 71N |43l. 1631] and 
AW [64] data (see, however, Ref. [65| for some criticism). However, for a more quantitative 
study of the constraints implied by pion production, and by nN- and iViV scattering data, 
a more complete and consistent analysis is necessary, which we postpone to a future work. 



V. SUMMARY AND OUTLOOK 



We performed a calculation of p-wave pion production amplitudes in A^V collisions in 
three different channels (pn —* pp7t~ , pp —* dix + and pp —> pnir + ) in the framework of chiral 
effective field theory. The relevant partial wave transition that depends on the (NN) 2 7r 
low-energy constant d is 3 S\ — > ^Sop for the first channel and ^So — > 3 S\p for the others. 
Therefore, it is clear that the study of different channels of the pion production reaction 
AW — > NNtt probes the corresponding operator in very different kinematical regimes and, 
thus, provides a non-trivial test for the validity of the employed approach. Our analysis of 
all the three channels resulted in values for the LEC d that are consistent with each other. 
In addition, we also point out an inconsistency in the partial wave analysis for pp — > pnii + 
carried out in Ref. [36|. Our findings can be interpreted as an indication that the source of 
the discrepancy reported in Ref. |3l|] is not due to the difference in the kinematics between 
AW — > NNn and tritium /3-decay, but rather caused by the inconsistency in the partial 
wave analysis for pp — > pri7c + as well as by some technicalities with respect to the work of 
Ref. jH] that we also discussed in our paper. 

Our investigation implies that calculations within effective field theory yield reliable re- 
sults for pion production in AW collisions utilizing the same value for d even though the 
corresponding contact term enters at very different kinematics in the reactions pn — > ppir~ , 
pp — ► dir + and pp — > pmr + . To confirm this conjecture, (i) one needs to reanalyse the re- 
action pp — > pnn + using the complete experimental information available for pp — > ppir as 
input and (ii) one needs new data for the process pn — * ppn~ at lower energies. We would 
like to stress that the near-threshold measurement of the reaction pn — > (pp)is 7i~ , where 
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(pp)i So signifies that the final pp state is constrained to be in the S wave by a kinematical 
cut, is the cleanest way to extract information on the contact term from pion production pro- 
cesses. Such measurements are already under way at COSY. Indeed, in the near future both 
pp —>■ {pp)is 7r ° an d pn —>■ {pp)^ 71 ^ wm be measured even with polarized initial state [39]. 
In addition, a consistent calculation for both tritium beta decay as well as low-energy pd 
scattering should be performed. We plan to perform these calculations in the future. 
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APPENDIX A: REACTION AMPLITUDES 

In this Appendix we present expressions for the matrix elements for the reactions we 



consider. To calculate them, we used the technique developed in Ref. [66 . 



1. General Considerations 

Let us consider pionic reactions involving the NN system, for example NN — > NNw, 
ird — > NN, etc. In the most general case, an amplitude corresponding to the matrix element 
of a particular production and/or absorption operator between states with given initial 
(j, I, s) and final (J 1 , 1', s') total angular momentum of a nucleon pair, its orbital momentum 
and total spin 2 is written as 

^l full [j7 S , j'Z's'] = ^ tree j'Z' S '] +^l FSI +^ ISI [j7 S , j'Z' S '] +^l ISI+FSI b7 S , j'Z's'] ,(A1) 

where "tree" stands for the tree production amplitude, i.e. where there is no NN (or NA) 
interaction both in the initial and in the final state, and FSI, ISI, ISI+FSI refer to the 
amplitudes with final state, initial state, and both final and initial state interaction included, 
in order. In this equation we imply that the spin-angular part (as well as the isospin part) 
of the amplitudes are factored out. Note that since there is a third particle that carries 
angular momentum, the pion, the total angular momentum j of the initial two-nucleon state 



2 In order to unambiguously specify the partial wave, the pion angular momentum should, in general, also 
be given. We, however, omit it since it is only the p-wave pion production that is considered here. 
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can be different from that of the final two-nucleon state, f . Obviously, the total angular 
momentum of the final particles has to be equal to that of the initial ones. Given the tree 
amplitude as a function of the initial p and final p' relative momenta, A tree [jls, j'l's'](p, p'), 
the remaining amplitudes are given by the following formulae: 

.FSir-, , n ,n _ v f d 3 q A^jjls, fl" s"}(p, q) M[j' , l"s", ls](q,p') 

,lSlr 7 ■/,//, V- f M[jJs,l"s"}(p,q)A^[jl"s",j'l's'}(q,p') 
A Uls,jls\-2^ I (27r)3 4M 1 M 2 [gV(2/i 12 )- J E-^0] ' lAdj 



J",s" 



7// _// 1111 ~fff J 



V's" l"'s'" ' 



d 3 q dH 
(2vr) 3 (2tt) 3 



M [j, Is, l"s"] (p, g) ^ [;7 V ; , jT ff ^ ff ] (g,£) M [f, l"'s>", l's'} (i, p') 
X 4M 1 M 2 [g 2 /(2/ii 2 ) -E-iO] •4Mi/M 2 /[£ 2 /(2/ii/ 2 /) -.B' -zO] ' 1 ' 

where M lj2 (Mi' i2 ') are the masses of the particles in the intermediate state that are related 
via the NN interaction to the initial (final) state, //i 2 {nvv) are the corresponding reduced 
masses, E (E') is the energy of the initial (final) two-nucleon state in its center-of-mass 
frame, JA\j,liSi,lfSf] is the iVTV half-offshell .M-matrix corresponding to a transition from 
the state (jhsi) to the state (jlfSf), and the sums are over all the intermediate states with 
given j, j', I, I', s, and s'. We use the following relation between the .M-matrix and the 
commonly used T-matrix: M. = — 8n 2 y/ MiM 2 M 3 M 4 T, where the Mj are the masses of 
interacting particles. 

The formulae given above also hold for the case when there is a transition through an 
intermediate iVA state going to a final (from an initial) state via an NN — NA interaction. 
In this case the NN .M-matrices have to be replaced by the appropriate NN — NA matrices, 
and also the propagators entering Eqs. (1A2I) - (1A4I) that correspond to the iVA intermediate 
state have to be modified according to 

(A5) 



AM 1 M 2 [q 2 /(2 f i 12 ) -E-iO] 4V2MiM 2 [? 2 /(2//i 2 ) - (E - AM) - iO] 

where AM is the nucleon-A mass difference (note also the factor l/y/2). Of course, a tree 
diagram with a iVA initial or final state gives a nonzero contribution only when it is inserted 
as a building block into those of FSI and ISI diagrams that have iVA as an intermediate 
state. 

In case of a deuteron in the final state, the corresponding M. matrices should be replaced 
by the deuteron wave functions according to 

*4 FSI b'^ 1] = -i^r E / 4% llV](p, g)/V"(g), (A6) 

V 2M N v , J [Any 

where ip 1 ' (q) are the deuteron wave functions corresponding to the angular momentum I", 
normalized by the condition 

^ ((^ (g)) 2 + (^ 2 (g)) 2 )=l. (A7) 



(27T) 
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Thus, the two-nucleon propagator for the deuteron in the final state is absorbed in the wave 
functions and the normalization has changed. Analogous expressions can be written down 
for the deuteron in the initial state and also for the deuteron in the initial and final states. 
Note that in the case of the deuteron in the inital and/or final state the tree diagrams appear 
only as building blocks for the calculation of the ISI/FSI and ISI+FSI diagrams according to 
Eqs. ( lA2l) - (jA4l) and IA6I) . respectively. They do not contribute independently because then 
there are no free nucleons in the initial and/or final state. 



2. The reaction pn — > ppir 

Here and below we use the spectroscopic notation 2S+l Lj for the NN and NA partial 
waves rather than the [jls] notation used in the previous section. The transitions that 
contribute to the reaction pn — > ppn~ at energies close to threshold are 3 Si — > x Sop, 3 D X — > 
l Sop in the isospin-zero initial state, and 3 Pq ~~ ► 1 SqS in the isospin-one initial state. The 
spin-angular structure of the amplitude reads 



M 



pn^ppir 



AxiSp) + d(SK) + C 2 S((pK)p - -K 



(A8) 



where S = X2 T^Xii — XvT^Xv denote normalized spin structures corresponding to the 

initial spin-triplet and final spin-singlet states, in order. Here and below, p,p', denote unit 
vectors of initial and final relative momenta of two nucleons and that of the pion momentum, 
respectively, and the with corresponding indices stand for the spinors of the initial and 
final nucleons. In turn, A±, C±, and C2 are the amplitudes corresponding to the 3 Po - > x Sqs, 
3 Si — > ^SqP, and 3 D 1 — >■ ^Sop transitions, in order. They are related to the corresponding 
amplitudes in the JLS basis via 



fullr3 



(A9) 

(A10) 
(All) 



The observables we consider are expressed in terms of these amplitudes in the following way: 



A,. 



d 2 a 



dVtdm 2 p 



d 2 a 
dVLdm 2 p 



Pcut 



4AM p 2 p (4vr) 4 M 7V sp 



K{p')p dp 



\M 2 + 



2C 

+2Re(A* l (C 1 + —l))cos9 w 



2Re (C*Ci 



Pcut 



4AM p 2 p (4vr) 4 M 7V sp 



K(p')p' dp 



\a 



9i 
3 

|2n 



cos 2 6 n 



(A12) 



C 2 



sin 26 n Im (C*C 2 ) - 2 sin 6 W Im {A\{C X - 



(A13) 



where p cut is the maximum relative momentum of the final protons in the measurements at 
TRIUMF 0, Bj, AM 2 p = (2M N +p 2 cut /M N ) 2 - (2M N ) 2 w 4p 2 cut , K(p') is the momentum 
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of the final pion, and s and p are the invariant energy squared and the relative momentum 
of the initial nucleons, in order. 

Below we give the expressions for the tree amplitudes «4[ 3 Si, ^So], A[ 3 D ly ^o] resulting from 
various pion production mechanisms as well as those for the relevant production amplitudes 
involving the A isobar. Note that from here on we suppress the label "tree" on the tree- level 
transition amplitudes. We also explain how we extract A\ from experimental data. 

a. Direct production 



AfS^S^p 1 ) 
AfD^So^p') 



C_ r dSl E 
,/2J 



-K + ^(P%) (2n) 3 5^(p'-p+k w /2) , 
M N 



- kMfKf- 1] + %r Wk) (pp') - (p%)} 

M N 



V2~J 4tt 

x(2n) 3 5^ 3 \p'-p+k n /2), 



where C 



. SM 2 N g A 

-l and uj„ 



(A14) 



(A15) 



k 2 + m 2 is the energy of the final pion. Here we included 



UV2 ' 

both the leading ttNN vertex and its recoil correction which enters at NNLO. 



b. Production via the A(1232) isobar 

The A(1232) contribution comes from the A^A intermediate states. In the reaction pn — > 
pp7r~ with the initial isospin of the pn system being 1 = 0, the NA <-> A^A" transitions are 
allowed only in the final state interaction. As we consider those kinematical configurations 
where the relative kinetic energy of the final protons is small, it is only the ^So final state 
that contributes. Therefore, the only coupled channel where the A(1232) contributes is 
5 Dq(NA) — > ^SoiNN). For p wave pions the relevant amplitudes that correspond to the 
3 Si(NN) -> 5 D (NA) and ^(NN) -> 5 D (NA) transitions in the production operator 
read: 

<hh [ - t[3(p'i) 2 -l] 



A[ 3 S^D Q ]{p,p') = C A J^ 
A[ 3 D^D,]{p,p') 



.8M N M A h A M N 
where C A = — i — — — — — i/— , and it 



V2~J 4tt 

x(27r) 3 5 (3 \p'-p+k^), 

M N 



the JVA state. 



3f n V2 |M A : 



(2n) 3 5^ 3 \p'-p+k^), (A16) 
K[9(p'p)(p%)(pk)-Z(ph) 2 -3(p%) 2 +l] 

(A17) 

and p' is the relative momentum of 



M N + M A 



c. Rescattering via the s-wave WT vertex 



AFS u %](p,i/) = 



C u^k w r dVL 



n f Uie'p 

~ J 177? 



(p-p') 2 + ml 



2(p — p 



l\2 



p 



l\2 



mi 



, (A18) 
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A\ 3 n 1 ?](nr i >) - fdnp 3(p'P-p) 2 -(p-p') 2 (A v 

5o](P ' P) - i^Tf i 17 [(p-p') 2 + "«] 2 • (M9) 

Note that in these expressions, and also in the expressions for the amplitudes that stem 
from operators with C3, C4, and recoil corrections to the WT vertex (see below), we keep only 
the leading term in the expansion in powers of k n /p. The same is true for the corresponding 
amplitudes in the reactions pp — > dn + and pp — > pnn + . 



d. Operators with C3, 04, and recoil corrections to the WT vertex 



AfD^S^p') 



AC K f dVLj;, 



3& {P-P' y2 



3/1 7 47T 
2C /■ cKlp> 



C 



p' 2 — p 2 



Z fly/2 J 4vr 



a 



(p — p') 2 + m2 8M N (p — p') 2 + m 2 

■tp ^jp'p-p) 2 - ip-p'f 
(p — p') 2 + ml 



,(A20) 



where C- Sl = ^ + 04 + -^ 



1 3((p / p) 2 -p 2 ) -p /2 +p 2 
4M^ (p-p') 2 + ml 

1 



C,-^ 1 — C3 — C4 



1.1/ v' 2 '" '" AM, 



N 



(A21) 



e. Contact term 



9C h 

A[%/S }(p, P >) = ±*Ld. 



A[ 3 D u %](p,p') = 0. 



9 a 



(A22) 
(A23) 



/. Coulomb interaction 

Since the final two protons are at low relative momenta, there are sizable effects from 
the Coulomb interaction between the two protons. The effect of the Coulomb interaction 
was taken into account along the lines of Refs. (67|, 68|. To be specific, we multiply all tree 
diagrams that do not contain the A isobar by the Gamow-Sommerfeld factor 



G{p') 
7(P0 



27T 7 (p') 



exp 27T7(j/) — 1 
M N 
lap 1 ' 



1/2 



(A24) 



where a is the electromagnetic fine structure constant. At the same time, the half-offshell pp 
M. matrix in the partial wave that we use in our calculation is corrected for the Coulomb 
interaction according to 



M cc {q,p') = M{q,p 



G(q) M^(p',p') 
G{p>) M{p',p') 



(A25) 
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where A4(q,p') and Ai cc (q,p') are, in order, the half-offshell pp M. matrices without and 
with the inclusion of the Coulomb interaction, whereas A4(p',p') and Ai cc (p' ,p') are the 



corresponding pp onshell Ai matrices. See Ref. [67| for more details. 

As far as diagrams with the A are concerned, where we have the transition 5 Dq(NA) — > 
^oiNN), we apply the following argument in order to take into account the Coulomb in- 
teraction. First, we note that the typical relative momenta of the intermediate A^A state 
are large so that the Coulomb interaction in this intermediate state is expected to be unim- 
portant. In order to take into account the Coulomb interaction between the protons in the 
final state, we multiply the amplitude with NA by the ratio of the inverse Jost functions 

F cc (p') 

GkaV) = W (A26) 

where F^ p c (p') and F pp {p') are the pp inverse Jost functions with and without Coulomb 
interaction, respectively. They are related to the corresponding M. matrices via 

F PP (P) - G(P) + ^J (27r)3 q2 _ p , 2 _ t0 , (A27) 

F ^ P ) = l + AM-J (2vr)3 g2 -^-,0- (A28) 

Notice that the overall normalization of the Jost functions is of no relevance, for they only 
occur in the ratio (however, both Jost functions have to have the same normalization factors). 



g. The contribution of the 1 = 1 initial state 



We write A 1 in Eq. (|A8[) in a form which takes into account the ISI phase as well as the 
dependence of the final state Jost function on the momenta: 

A 1 = iXexpti^F^V). (A29) 

The constant (real) factor X is adjusted in such a way to reproduce the corresponding 
amplitude extracted from the TRIUMF data using a partial wave analysis 0, 47]. The sign 
of X is adjusted to the behaviour of observables in pn —>■ ppiT~ . 



3. The reaction pp — ► d7r+ 



The transitions that contribute to the reaction pp — ► dii + at energies close to threshold 



are % - 3 S lP , % - 3 D lP , l D 2 - 3 S lP , X D 2 
spin-angular structure of the amplitude reads 



3 D lP , and 3 Pi -> 3 Sis, 3 Pi -> 3 D lS . The 



M 



pp— >dn+ 



C (S xp)e+ dl{ke) + C 2 l[{pk){pe 



- 3 m\ 



(A30) 



where e is the deuteron polarization vector, S = xI^&Xii % = xI^Xi are normalised 
spin structures corresponding to the initial spin-triplet and spin-singlet states, in order. 
Here, the x's refer to spinors of the initial nucleons, and C , C\, and C 2 are the amplitudes 
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corresponding to the 3 P\ — > es, — > ep, and 1 D 2 ^ ep transitions, in order. Further, we 
denote the final states as el, where e stands for the deuteron final state, and I is the angular 
momentum of the final pion relative to the deuteron. The amplitudes Co, C\, and C 2 are 
related to the corresponding amplitudes in JLS basis via 



Co 
C 2 



^ fuU [ 3 Pi,l], 



^15 



^2,1]- 



(A31) 
(A32) 

(A33) 



Note that, as it can be seen from Eq. (]A6|) . the amplitudes of the transitions to the deuteron 
state are sums of the amplitudes where the transition goes to the S'-wave component and 
those where the transition goes to the D-wave component of the deuteron wave function. 
Note also that here the total amplitudes are distinguished only by the initial state, as 
the final state, the deuteron, is the same for all transitions. However, at the level of tree 
amplitudes, one has to distinguish between the ones that correspond to transitions to the 
.S'-wave component and those to the .D-wave component. 

The observables under consideration are expressed through the amplitudes C , C\, and 
C2 as 



.4 



da 

dn 

da 



2567T 2 SJ9 



2|Cnl 2 



Ci| 2 + i|C 2 | 2 (3cos 2 W + 1) + ^Re(CiC*)(3cos 2 6 n - 1) 
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2 sin 6 n cos Im (Cq (Ci 



Co 



))■ 



(A34) 
(A35) 



y dn ~ 256vr 2 sp r — ^-^-^ 3 

Below, we give expressions for the tree amplitudes A[ 1 D 2 , 3 Di], A[ 1 D 2 , 3 Si], A^Sq^Di], 
A^So, 3 Si), contributing to C\ and C 2 , as well as those for the relevant production amplitudes 
involving the A isobar. We also provide details of the determination of C . 



a. Direct production 



A^SA^p') 
AfSoM^p') 

a^dA}^') 



c 



dQ.i 



47T 

C r 



-K+^r (P%) (2^) 3 5 (3) {f -p+ K/2) 
M N 



V2J 4vr 
x(27t) 3 5 (3 \p'-p+K/2) 



K[3(p'K) 2 -l} + ^(p%) 

M N 



dflj 



10 J 4tt 

x(2n) 3 5 {3 \p'-p+K/2) 
C f3~ rdfl 



K[z(pK) 2 -i] + 



M, 



N 



'){pK)-p'K\ 



^7T 



(A36) 



(A37) 



(A38) 



-K[%p'K){p'p){pK)-3{pKf-3{p'K) 2 +i} 



V2V10 7 4tt 

+^[3(pp')(ih)-p'h} {2^) 3 5^\p'-p+K/2) . 



M, 



A' 



(A39) 
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b. Production via the A(1232) isobar 



In the reaction pp — ► dii + , the initial isospin of the pp system is / = 1, so the 
NA «-> AW intermediate states are allowed only in the initial state interaction. The 
coupled channels that contribute are 3 %{NN) -> 5 D (NA), ^(NN) -> 5 D 2 (NA), 
and 1 D2(NN) — > ^(iVA). The relevant transitions in the production operator for p- 
wave pions are 5 D (NA) -> 3 Si(NN), 5 D (NA) -> 3 Dx(NN), 5 D 2 (NA) -> ^(NN), 
5 D 2 (NA) -> ^(N^fS^NA) -> 3 Si(NN), and 5 S 2 (NA) -> 3 D X (NN). The expressions 
for the corresponding amplitudes read: 



^[ 5 d , 3 ^i](p,pO = 



/i'77 



47T 



M3(p%r) 2 - 1] 



^[ 5 iV£i](p,p') = c- A W- 



v/2 7 4tt 

x(27r) 3 ^(p"-p+4/2) 

2T /■ (in 



20 J 4vr 
6" /• c% 



(2tt) 3 (5^(p , -P+^/2), (A40) 
fc* [9(p p) (p 'Mr) (p4) - 3 (p^) 2 - 3(p '^) 2 + 1] 

(A41) 

(2vr) 3 (5^(p , -p+4/2), (A42) 



-Av[3(p^) 2 -1] 
-K[9(z>pf-6(pp')(z>kWK) 

1 11, 



+2(p^) 2 + 2(p'£;, 



(27r) 3 ^(p"-p+^/2), 



(A43) 
(A44) 



^[ 5 S 2 , 3 £i](p,p') 



"10W 17 



^[3(j5'i) 2 -l] (27r) 3 ^(p / -p+4/2). (A45) 



c. Rescattering via the s-wave WT vertex 



AfSo^&p') 
A[ x D 2 3 S x ]{p,p') 

ai'd^d.Up') 



C uj n k n r dQ 



2/2 J 4vr (p-pif + ml 



1 - 



2(p — p 



l\2 



p') 2 + m 2 



p') ~~ (p ~~ p 



A2 



3/2^/2 7 4tt [(p-p') 2 + m 2 ] 2 
C^K fdQ f 3(p'p-p) 2 -(p-p') 2 



tfVSOJ 47T 

3 C Ld-n-k-K f dVtfi 



— p') 2 + m 2 
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l\2 



V2 f 2 V30 J 4tt (p-p') 2 + m 2 
9(p—p'p) (pp' — p' ) (pp') — 3 (p — p 'p) 2 — 3 (p' — pp' ) 2 + (p — p 
3[(p — p') 2 + m 2 ] 



(A46) 
(A47) 
(A48) 



/\2 



.(A49) 



3 We do not take into account the channel 1 D2{NN) — > 3 Z?2(-/VA) because the corresponding AW — > TV A 
.M-matrix is subleading according to the power counting and it is also numerically small at the energies 



considered 
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d. Operators with c%, 04, and recoil corrections to the WT vertex 



A[%,%}(p,p') 



A^D^D^p') 



ACK r da 



2Ck 

71 



All 

dQ f 



a 



P 2 — p' 2 



1 3((pp') 2 



C 



3 D, 



(p — p') 2 + m 2 8M N [p — p') 2 + m 2 
3{pp' -p') 2 - {p-ff 



,(A50) 



(p — p') 2 + ml 



P 



f2\ 



p 2 + p' 2 



AM, 



N 



(p — p') 2 + m 2 



c 



3 Di Hp'p-p) - ip-p') 



A2 



(p — p') 2 + m 2 . 



1 3((p'p) 2 -p 2 ) -p' 2 +p 



AM 



N 



(p — p') 2 + ml 



f 2 VTEJ 47T 



dQf 



C 



3 Dl j9(p-p'p)(pp'-p'){pp') 



(p — p') 2 + ml 



-3{jp—p'p) 2 — 3{p'—pp' ) 2 + {p—p 
(p — p') 2 + ml 
1 [9 (p'p +p)(p'— pp' ) (pp' ) 



i\2 



AM N \ {p-p') 2 + ml 

—3(p'+pp')(p'—pp')—3(p'p + p)(p'p — p) + p' 2 —p 2 
(p — p') 2 + ml 

1 \ [3(pp , ) 2 -l](p-p / ) 2 " 
AM N ) (p-p') 2 + ml 



+3 G4 + 



(A51) 



(A52) 



(A53) 



e. Contact term 



9C h 

A[ l SA}{p,p') = ^^d, (A54) 

9a 

A[%M]M) = 0, (A55) 
AfD^fap') = 0, (A56) 
A^D^D^p') = 0. (A57) 

/. Coulomb interaction 

For this reaction, the Coulomb interaction between the initial protons gives only a small 
effect since the initial relative momentum is large. On the contrary, the deuteron and pion 
in the final state are at low relative momenta, and the Coulomb interaction between them is 
taken into account at the level of experimental data, by factoring out the Gamow-Sommerfeld 
factors that stem from the ir + d Coulomb interaction 
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g. Pion production in the s-wave 



In order to calculate the amplitude of the s-wave pion production, that is, the value of 
Co, we took the result of Refs. (l7|, l69|, where the amplitude of the s-wave pion production 
was calculated up to NLO. We correct this amplitude by a constant factor so as to get at 
threshold the value of s-wave production parameter a = 252/ib. This quantity, a, is the 
total cross section divided by the final pion momentum in units of the pion mass, and the 
given value is extracted from the width of pionic deuterium [50|, U . 



4. The reaction pp — > pnir 



+ 



The transitions that contribute to the reaction pp — > pnir + at energies close to threshold 
are % -> 3 S lP , % -> 3 D lP , X D 2 -> 3 S lP , l D 2 -> 3 D lP , and 3 P X -> 3 Sis, 3 P X -> 3 £> lS . At 
very low energies, the final 3 .Di state does not contribute, however, the transitions to the 3 Di 
state contribute via the 3 Di «-> 3 Si coupled channel. At these low energies, the spin-angular 
structure of the amplitude reads 



M 



pp— ►pmr -1 



C (Sxp)S' + C{I{kS') + C 2 l[{fk) ipS') - ~{kS') 



(A58) 



where S = xi^&Xii 1 — xi^kXii = Xv&t^Xz' are normalized spin structures corre- 
sponding to the initial spin-triplet, initial spin-singlet, and final spin-triplet states, in order. 
Here, Co, Ci, and C 2 are the amplitudes corresponding to the 3 P\ — > 3 S\s, ^So — > 3 S\p, and 
1 D 2 — > 3 5ip transitions, in order. Their relation with the corresponding amplitudes in the 
JLS basis is given by 



C x = A m ?So,% 



C 2 = \jf A m ^%]- 



(A59) 
(A60) 
(A61) 



Besides these amplitudes, there are also amplitudes that correspond to the transition to the 
isosp in-one ^So final pn state. However, these amplitudes do not interfere with C , C\, and 
C 2 because of different spins in the final state and generate only an additive part in the cross 
section — see also the discussion in the text. The observables in the reaction pp — ► pmr + 
can be expressed through C , C\, and C 2 in the following way: 



da 



da 



dQ 2(4vr) 4 M 7V sp 



k n (p')p' 2 dp' 



2|C | 2 + |C 1 | 2 + i|C 2 | 2 (3cos 2 ^ + l) 
y 



+-Re (CiC 2 *)(3cos 2 6„-l) 



+ 



da 

In 



1=1 



Pmax 



A ■ 

v dVL 2{4^) 4 M N sp 



k n {j)')p dp ■ 2 sin 0^ cos </> Im(Co(Ci 



(A62) 
(A63) 



24 




FIG. 9: Leading diagrams for ttN scattering with intermediate Delta. Shown are the s-channel (a) 
and the u-channel (b) contribution. 

da 1=1 

Here, p max is the maximum relative momentum of the final nucleons, and -— is the 

dll 

contribution of the final pn state with isospin one to the cross section. 

The expressions for the transition matrix elements relevant for the calculation of C\ and 
C*2 are the same as for the reaction pp — > dir + , and therefore we refer the reader to the 
corresponding formulae, given in Sec. IA 31 above. However, we should make a remark about 
the value of Cq. Instead of calculating it at NLO we extracted Cq directly from the data - 
see the discussion in the main text. 



APPENDIX B: RESONANCE SATURATION VS. EXPLICIT DELTA 

In the chiral limit the masses of both the A and the nucleon stay finite and differ from each 
other. Thus, there is a well defined limit of QCD where m n /AM, with AM = — M^, is 
a small parameter. Consequently, the Delta degrees of freedom may be integrated out and 
the effects of the A isobar are then absorbed into the LECs of the resulting Lagrangian. 
On the other hand, when the Delta degrees of freedom are included dynamically in the NN 
system, certain selection rules apply. In particular, the A is allowed to contribute only, if 
the NN system has the total isospin equal one. It is instructive to discuss the implications 
of these selection rules for the reaction iViV — > NNtt using both the EFT with and without 
explicit Delta degrees of freedom. The corresponding discussion for ttN scattering can be 
found in Ref. (ilf f 

For simplicity, let us start from elastic ttN scattering. Using the interactions defined in 
the main text one finds straightforwardly for diagram (a) of Fig. [9] in the limit AM — > oo 



M N h\ 



.18/2AM 



q'tfjiZSij - ie ijk a k ){25 ba - ie bac r c ) , (Bl) 



where a (b) and q (q') denote the isospin quantum number and momentum of the incoming 
(outgoing) pion. Analogously, one finds for the u-channel diagram 

( M h 2 \ 

= ~i ( 18 p A A M J ViH( 2S ij ~ ie jik a k )(25 ba - ie abc r c ) . (B2) 

In the second line of the above equation we used Eqs. Thus, diagrams (a) and (b) 
are individually given by four terms. However, two of them get canceled when adding the 
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contributions together. The remaining two terms in the sum exactly resemble the structure 
of the C3 and C4 terms of the effective Lagrangian. One then finds for the Delta contribution 



to the LECs Cj the following result 41 



-2d 



h 2 
9AM 



(B3) 



Obviously this matching is only possible after adding the two diagrams together. On the 
other hand, the above mentioned selection rules for iViV — ■> NNtt are operative at the level 
of the individual diagrams. We now show how these two facts can be realized simultaneously. 

In full analogy to the expressions given above, one finds for the amplitudes corresponding 
to the diagrams of Fig. [2] involving the Delta: 



%A 



(c) 



' M N9Ah\ \ (2) / (2) s 1 / 

18f*AMJ a ^ q K 2 -ml qA \ ' 

' M 2 N g A h\ s 
18/1 AM, 



(i) 



(26; 



ba ^^-bacT, 



(1) 



2<l'i<l3 

q* — m* 



X 



%A 



(j.) 



(M 2 N g A h\ \ _ (2) / .,, 



\ 18/3 AM 

/ M|g^ \ (2) / (2) _ x 1 

I 18/3 AM I a V° y y«2_ 



2z {ey fc cr^ 1) 5 ab + 5ije 6ac r c (1) |) 



(B4) 



6a l^-abcT, 



(1) 



q* — m 2 



x 



2i [e ijk a^5 ab + f^e^cT^}) , 



(B5) 



where in both amplitudes the momentum of the outgoing (virtual) pion is labelled as q' (q). 
In both amplitudes the terms in [...] exhibit the same spin-isospin structure as the Cj terms 
discussed above, while those in {...} have a different structure. 

Therefore, when diagram (c) and (d) are added together, only the structures of the q 
parameters survive. However, due to the above mentioned selection rules, the diagrams 
contribute individually to different channels of NN —>■ NNtt. To be specific, the diagram 
(d) vanishes for the isospin-one to isospin-zero transition, i. e. for pp — > (d/pnj = o)ir + , whereas 
the diagram (c) does not. Furthermore, once the isospin matrix element is evaluated for 
the diagram (c) it turns out that the expression in {...} gives the same contribution as the 
one from [...]. The same holds for the isospin-zero to the isospin-one transition, i. e. for 
pn — > ppiT~ - the diagram (c) does not contribute whereas the terms in brackets {...} and 
[...] for the diagram (d) are equal. Therefore, in the limit AM — > 00, one indeed observes 
both properties simultaneously, namely that the iVA intermediate state does not contribute 
if the external NN state is in the isospin-zero state and that the Delta effects can be absorbed 
in local counter terms, namely C3 and C4. 

We are now also in the position to see how the pattern changes when we start to move away 
from the limit AM — > 00. Then the factors 1/AM that appear in Eqs. (1B1I) - (1B5I) should 
be replaced by the dynamical iVA propagators. One finds for the resulting combination of 
the propagators for both ttN scattering and NN — > NNtt 



1 



AM 



1 

AM + cj„ 



(B6) 
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where we used that ~ E tot and dropped terms significantly smaller than m^. In this 
expression the upper (lower) sign refers to the combination of propagators relevant for the 
terms that can (cannot) be mapped onto the Cj. Thus, the additional terms are suppressed 
by 

5 = u n /AM . (B7) 

Numerically S is already as large as 0.5 at threshold and grows as one goes to higher energies. 
Clearly, near the two-pion production threshold, 5 ~ 1 and it is necessary to keep the Delta 
as dynamical degrees of freedom. See Ref. |70J for a power counting that allows one to also 
study the latter regime. 
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